Quasiclassical fluctuations of the superconductor proximity gap in a chaotic system 
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We calculate the sample-to-sample fluctuations in the excitation gap of a chaotic dynamical system 
coupled by a narrow lead to a superconductor. Quantum fluctuations on the order of magnitude of 
the level spacing, predicted by random-matrix theory, apply if te <§C h/ Et (with te the Ehrenfest 
time and Et the Thouless energy). For te ^ H/Et the fluctuations are much greater than the level 
spacing. We demonstrate the quasiclassical nature of the gap fluctuations in the large-r^ regime by 
correlating them to an integral over the classical dwell-time distribution. 
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The universality of statistical fluctuations is one of the 
most profound manifestations of quantum mechanics in 
mesoscopic systems 1] . Classically, the conductance g of 
a disordered metal (measured in the fundamental unit 
2e 2 /h) would fluctuate from sample to sample by an 
amount of order (l/L) 3 / 2 <C 1, with I the mean free path 
and L the length of the conductor 0. Quantum me- 
chanical interference increases the fluctuations to order 
unity, independent of disorder or sample length. This 
is the phenomenon of universal conductance fluctuations 
0, Q . The same universality applies to a variety of other 
properties of disordered metals and superconductors, and 
random-matrix theory (RMT) provides a unified descrip- 
tion |5j . 

Chaotic systems (for example, a quantum dot in the 
shape of a stadium) share much of the phenomenology of 
disordered systems: The same universality of sample-to- 
sample fluctuations exists 0, Q, @ ■ What is different is 
the appearance of a new time scale, below which RMT 
breaks down 0, ^j. This time scale is the Ehrenfest 
time te, which measures how long it takes for a wave 
packet of minimal size to expand over the entire avail- 
able phase space. If te is larger than the mean dwell time 
Try in the system (the reciprocal of the Thouless energy 
Et = K/2td), then interference effects are inoperative. 
A chaotic system with conductance g x 2e 2 /h, level spac- 
ing 6, and Lyapunov exponent A has rrj = 2nh/gS and 
te = A -1 In (gro/ro), with tq the time of flight across the 
system pj|. The defining characteristic of the Ehrenfest 
time is that it scales logarithmically with h, or equiv- 
alently, logarithmically with the system size over Fermi 
wavelength [l2^ . 

The purpose of this paper is to investigate what hap- 
pens to mesoscopic fluctuations if the Ehrenfest time be- 
comes comparable to, or larger than, the dwell time, so 
one enters a quasiclassical regime where RMT no longer 
holds. This quasiclassical regime has not yet been ex- 
plored experimentally. The difficulty is that te increases 
so slowly with system size that the averaging effects of 



inelastic scattering take over before the effect of a finite 
Ehrenfest time can be seen. In a computer simulation in- 
elastic scattering can be excluded from the model by con- 
struction, so this seems a promising alternative to inves- 
tigate the crossover from universal quantum fluctuations 
to nonuniversal quasiclassical fluctuations. Contrary to 
what one would expect from the disordered metal 0], 
where quasiclassical fluctuations are much smaller than 
the quantum value, we find that the breakdown of uni- 
versality in the chaotic system is associated with an en- 
hancement of the sample-to-sample fluctuations. 

The quantity on which we choose to focus is the exci- 
tation gap £o of a chaotic system which is weakly cou- 
pled to a superconductor. We have two reasons for this 
choice: Firstly, there exists a model (the Andreev kicked 
rotator) which permits a computer simulation for sys- 
tems large enough that te ^ Try- So far, such simula- 
tions, have confirmed the microscopic theory of Ref. [ll| 
for the average gap (so) 0. Secondly, there exists a 
quasiclassical theory for the effect of a finite Ehrenfest 
time on the excitation gap and its fluctuations [l4| . This 
allows us to achieve both a numerical and an analytical 
understanding of the mesoscopic fluctuations when RMT 
breaks down. 

We summarize what is known from RMT for the 
gap fluctuations 0. In RMT the gap distribution 
P(eq) is a universal function of the rescaled energy 
(eo — Eg)/A g , where E g — 0.6 Et is the mean-field en- 
ergy gap and A g = 0.068 g 1 ^ 3 S determines the mean level 
spacing just above the gap. The distribution function 
has mean (eq) 
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1.21 A g and standard deviation 
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The RMT predictions for P(sq), in the regime te <§; 
tu, were confirmed numerically in Ref. |l3| using the 
Andreev kicked rotator. 

We will use the same model, this time focusing on 
the gap fluctuations Seq in the regime te ^ Tr>. The 
Andreev kicked rotator provides a stroboscopic descrip- 
tion (period To) of the dynamics in a normal region of 
phase space (area Mh c s) coupled to a superconductor in 
a much smaller region (area Nh c ff, 1 <C N <C M). We re- 
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FIG. 1: Root-mean-square value Seo of the gap divided by the 
RMT prediction fertMT , as a function of the system size M for 
dwell time M/N — 5 and kicking strength K = 14. The data 
points result from the numerical simulation of the Andreev 
kicked rotator. The solid line has slope 2/3, indicating that 
Seo depends only on M/N and not on M or N separately in 
the large-M regime. 



fer to this coupling as a "lead" . The effective Planck con- 
stant is h c g — 1/M. The mean dwell time in the normal 
region (before entering the lead) is td = M/N and the 
corresponding Thoulcss energy is Et — N/2M. We have 
set To and h equal to 1. The dimensionless conductance of 
the lead is g = N. The product 5 = AnEr/g = 2-k/M is 
the mean spacing of the quasi-energies e m of the normal 
region without the coupling to the superconductor. The 
phase factors e %Zm (m = 1,2,..,M) are the eigenvalues 
of the Floquet operator F, which is the unitary matrix 
that describes the dynamics in the normal region. In the 
model of the kicked rotator the matrix elements of F in 
momentum representation are given by [l6| 



Q n 



-(iir/2M)(n 2 +m 2 ) 



M -l/2 e (2iri/M)nm 



-(iMK/2-rr) cos (2-nn/M) 



(2a) 
(2b) 
(2c) 



The coupling to the superconductor doubles the di- 
mension of the Floquet operator, to accomodate both 
electron and hole dynamics. The scattering from elec- 
tron to hole, known as Andreev reflection, is described 
by the matrix 



1 - (1 - \V2)P T P 



pl/2 : 

with the projection operator 



-i±V2P T P 



-i\V2P T P 
1 - (1 - \V2)P T P 



(P T P) = s n 

v / rim 



1 if L<n<L + N 
otherwise. 



,(3) 



(4) 



Since we work in momentum representation, the lead de- 
fined by Eq. Q is a strip in phase space of width N 
parallel to the coordinate axis. One could alternatively 
consider a lead parallel to the momentum axis, if one 
would work in coordinate representation. We do not ex- 
pect any significant differences between the two alterna- 
tives. Putting all this together we arrive at the Floquet 
operator of the Andreev kicked rotator , 



T = V 1/2 



F 
F* 



v i/2_ 



(5) 



The matrix T can be diagonalized efficiently using the 
Lanczos technique in combination with the Fast-Fourier- 
Transform algorithm [17|. This makes it possible to cal- 
culate the quasi-energies e m and eigenfunctions ip m for 
systems of sizes up to M = 5 • 10 5 . The gap value eq is 
given by the eigenphase of T closest to zero. 

The Floquet operator JSJ) provides a stroboscopic de- 
scription of the electron and hole dynamics, which is be- 
lieved to be equivalent to the true Hamiltonian dynamics 
on long time scales t 3> tq. The support for this comes 
from two sides: (i) In the absence of superconductivity, 
and for varying parameters K and ?l e gf, the 1-D kicked 
rotator correctly reproduces properties of localized [Isj . 
diffusive [ljj , and even ballistic [2(j quasiparticles in dis- 
ordered media, (ii) In the presence of superconductivity, 
the kicked Andreev rotator, and extensions thereof, ade- 
quately describe quantum dots in contact with a super- 
conductor [13|, and give a proper description of quasi- 
particles in dirty d-wave superconductors |2l|. Since we 
will be giving a classical interpretation of our results, we 
also describe the classical map corresponding to the An- 
dreev kicked rotator. The map relates the dimensionless 
coordinate x n £ (0, 1) and momentum p n 6 (0, 1) at time 
(n + 1)tq to the values at time utq: 

Pn+i = p n ±(K/2Tr)sm[2ir(x n ± —)], (6a) 



Pn Pn+l 

2 2 



(6b) 



The upper and lower sign correspond to electron and 
hole dynamics, respectively. Periodic boundary condi- 
tions hold both for x and p. The quasiparticle reaches the 
superconductor if \p n +\ — Pioadl < N/2M, where pi ca d is 
the center of the lead. At the next iteration the electron 
is converted into a hole and vice versa. 

We study a system with kicking strength K = 14 (fully 
chaotic, Lyapunov exponent A = 1.95) and vary the level 
spacing 6 = 2-k/M at fixed dwell time T£> = M/N = 5. 
Sample-to-sample fluctuations are generated by varying 
the position pi ea d of the lead over some 400 locations. 
The resulting M-dependence of Seq is plotted in Fig. 
on a double logarithmic scale. We have divided the value 
5eq resulting from the simulation by the RMT prediction 
feRMT from Eq. QJ. The numerical data follows this 
prediction for M < 10 3 , but for larger M the fluctuations 
are bigger than predicted by RMT. For M > 10 4 the ratio 
5eq I femiT grows as M 2 / 3 (solid line). Since #£rmt oc 
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FIG. 2: Left panels: Husimi function J7J for the electron component of the ground-state wavefunction tpo of the Andreev kicked 
rotator, for two different positions of the lead. The parameters are M = 131072, td = M/N = 5, K = 14. The calculated 
values are scaled by a factor 0.019 (0.017) in the top (bottom) panel, so that they cover the range (0, 1), indicated by the gray 
scale at the top. Right panels: The corresponding classical density plots of all trajectories which have a time t > 7 between 
Andreev reflections. The calculated values are rescaled by a factor 0.30 (0.32) in the top (bottom) panel. 



M~ 2 / 3 , this means that Seq is independent of the level 
spacing S — 2tt/M at fixed dwell time td — M/N . This 
suggests a quasiclassical explanation. 

To relate the fluctuations of £o to the classical dynam- 
ics, we first examine the corresponding wavefunction i/jq. 
In the RMT regime the wavefunctions are random and 
show no features of the classical trajectories. In the qua- 
siclassical regime te > td we expect to see some classical 
features. Phase space portraits of the electron compo- 
nents of the wavefunctions are given by the Husimi 
function 

H(n x ,n p ) = \{il> e m \n x ,n p )\ 2 . (7) 

The state \n x ,n p ) is a Gaussian wave packet centered at 
x = n x /M, p — rip/M. In momentum representation it 
reads 

(n\n x , Up) OC e -An-n p f/M e 2^n/M ^ (g ) 



In Fig. |21 left panels, the Husimi function of i/jq is 
shown for two lead positions. Shown is a logarithmic 
gray scale density plot of the Husimi function, with light 
(dark) areas corresponding to low (high) density. The 
lead is visible as a light strip parallel to the x-axis. It is 
clear that these wavefunctions are not random. We ex- 
pect that the structure that one sees corresponds to long 
classical trajectories, since the wavefunctions are for the 
lowest quasi-energy. To test this expectation, we show in 
the right panels (on a linear gray scale) the correspond- 
ing classical density plots for all trajectories with dwell 
time t > t* . A total of 3 ■ 10 5 initial conditions (xo,Po) 
for these trajectories are chosen uniformly in the lead. 
Each new iteration of the map JSJl gives a point (x n ,p n ) 
in phase space, which is kept if the time of return to the 
lead is greater than t*. We take t* = 7, somewhat larger 
than the Ehrenfest time t e = A -1 In (N 2 /M) = 4.4. The 
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FIG. 3: The data points (left axis) are the quantum me- 
chanical gap values eo of the Andreev kicked rotator as a 
function of the position pi ca d of the lead, for parameter val- 
ues M = 131072, t d = M/N = 5, K = 14. The solid 
line (right axis) is the reciprocal of the mean dwell time 
(t}« = CT tP[t)dt/ /".„, P(t)dt of classical trajectories longer 
than t* = 7. 

plot is not particularly sensitive to the value of t* , as long 
as t* > te- There is a clear correspondence between the 
quantum mechanical Husimi function and the classical 
density plot. We conclude that the wavefunction of the 
lowest excitation covers predominantly that part of phase 
space where the longest dwell times occur. 

To make this more quantitative we show in Fig. [21 the 
gap value from the quantum simulations as a function of 
the lead position. The solid curve results from a classical 
calculation of the mean dwell time of those trajectories 
with t > t*, for the same value t* = 7 used in Fig. |21 



More precisely, it is a plot of 

(t). S?tP(t)dt' [> 

with P(t) the classical dwell time distribution. We see 
that the sample-to-sample fluctuations in the gap Eq cor- 
relate very well with the fluctuations in the sample-to- 
sample mean dwell time of long trajectories. Again, the 
correlation is not sensitive to the choice t* > Te- Such a 
correlation is in accord with recent theoretical work |14| , 
in which an effective RMT description is expected to hold 
for the part of phase space with dwell times greater than 
the Ehrenfest time. But we should emphasize that the 
agreement is only qualitative. In particular, the relation 
£o ~ l-5/(i)* — 0.07 that we infer from Fig. 3 is different 
from the relation e = 0.3/(i)* that would be expected 
from RMT. While the theory of Ref. H2 has been found 
to be in good agreement with the average gap value (eo) 
[l3| . it is not clear how it compares to the data of Fig. 3. 

In conclusion, we have investigated the transition from 
quantum mechanical to quasiclassical gap fluctuations in 
the superconductor proximity effect. The transition is ac- 
companied by a loss of universality and a substantial en- 
hancement of the fluctuations. Our numerical data pro- 
vides qualitative support for an effective random-matrix 
theory in a reduced part of phase space [l4|, as is wit- 
nessed by the precise correlation which we have found 
between the value of the gap and the dwell time of long 
classical trajectories (see Fig. 3). It would be of interest 
to investigate to what extent quasiclassical fluctuations 
of the conductance in a ballistic chaotic system are sim- 
ilar or different from those of the superconducting gap 
studied here. 
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